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£ • Abstract 
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0^ ■ We obtain bounds on the decay exponent A of the autocorrelation function 

; 

in phase ordering dynamics (defined by liim^^ ( (j)(r, tx)(j)(r, t^)) ~ L(t,2)~ x ) 
. For non-conserved order parameter, we recover the Fisher and Huse inequality, 

A > d/2. If the order parameter is conserved we also find A > d/2 if t\ = 0. 
However, for t\ in the scaling regime, we obtain A > d/2 + 2 for d > 2 and 
A > 3/2 for d = 1. For the one-dimensional scalar case, this, in conjunction 
with previous results, implies that A is different for t% = and t% ^> 1. In 
2-dimensions, our extensive numerical simulations for a conserved scalar order 
parameter show that A ~ 3 for ti = and A ~ 4 for t\ S> 1. These results 
contradict a recent conjecture that conservation of order parameter requires 
A = d. Quenches to and from the critical point are also discussed. 

PACS: 64.60.My, 68.35.Fx 
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Phase separation dynamics proceeds when a system is quenched from its high tem- 
perature, homogeneous phase to a low temperature, inhomogeneous phase (where several 
phases coexist in equilibrium). Due to its simple description yet rich behavior, phase order- 
ing dynamics has greatly enhanced our understanding of non-equilibrium processes ]Ij|. At 
late times, the spatial distribution of domains can be described by a single time- dependent 
length, L{t) which typically grows algebraically in time, L(t) ~ t 1 ^. This results in a scale 
invariant equal-time correlation function C(r,t). More recently it has been realized that 
the unequal-time correlation function is scale covariant. In particular, the asymptotic decay 
of the two-time autocorrelation function, C(r, ti,^) = ( <K r ) £i)0(O, £2)) defines an indepen- 
dent exponent A, via lim tl<i2 C(0, t 1 , t 2 ) ~ (L(ti) / L(t 2 )) x ■ This exponent bears no apparent 
relation to the growth exponent z and so its value provides a sensitive test for approximate 
theories of phase ordering kinetics Although the autocorrelation function has been 

studied extensively for non-conserved order parameter dynamics [ppHiO|, there has been 
hardly any work on conserved dynamics. However in a recent Letter, Majumdar et al. have 
shown numerically and analytically that that A = 1 for m — 1, d = 1 and ti = where m is 
the number of components in the order parameter TI]. It has been further argued JTTIJT^ 



that the conservation of order parameter demands that A = d for all m. 

In this Letter, we obtain lower bounds on the decay exponent A. For non-conserved 
order parameters, A > d/2 independent of £1, consistent with a general argument of Fisher 
and Huse 0. For conserved order parameters, we also obtain A > d/2 for ti = (assuming 
the quench is from a high temperature phase). However, for ti in the scaling regime, we 
find that A > d/2 + 2 for d > 2 and A > 3/2 for d = 1. This difference arises from the 
small k behavior of the scattering intensity S(k,ti). In conjunction with the exact result 



for A, for the 1-dimensional scalar model, with ti = [|TTJ, we conclude that for d — 1, A 
depends on whether £1 = or t\ ^> 1. To carry out the investigation in higher dimensions, 
we perform an extensive numerical integration of the Cahn-Hilliard equation (see Eq. (|J) 
below) in d = 2. We find that A ~ 3 for t\ = and A ~ 4 for t\ in the scaling regime. This 



is inconsistent with the recent conjecture that A = d |TT|JTM1- We discuss why this conjecture 



fails. We also derive bounds on A for quenches to and from the critical point. Our results 
easily extend to vector order parameters. 

We begin by obtaining the lower bounds on A. The equal point auto-correlation 
C(ti,t 2 ) = C(0, ti,t 2 ) is related to the k space auto-correlation S(k,ti,t 2 ) by 

C(t u t 2 ) = Jdk( 6Mti) ty-kftj) }= Jdk S(k, £1, h). 

Here <p(r, t) is the order parameter at point r and time t and 5<p(r, t) = <p(r, t) — m with 
mo = V" 1 J dr (p(v,t) and the fourier transform <5</> k (£) = V^ 1 ^ 2 J dr e _4kr 50(r, t). The 
angular brackets indicate an average over initial conditions. For a critical quench ( mo) = 
and (itLq) is O^V^ 1 ), whereas for an off-critical quench (m ) is 0(1). 
Using the Cauchy-Schwartz inequality, we find 

C(h,t 2 )< J dk(50 k (£ 1 )^_ k (£ 1 ))V 2 (^ k (£ 2 )50_ k (£ 2 )) 1 / 2 , 

~ J dkSfctjWSfch) 1 ' 2 . (1) 

where S(k,t) = S(k,t,t). 

Now assume t 2 to be in the scaling regime with t 2 > t\. At late times, the scattering is due 
to the sharp interfaces or defects. The k modes, <50 k at times t\ and t 2 will be uncorrelated 
when the interfaces move a distance greater than 27r/fc so that S(k,ti,t 2 ) decreases rapidly 
for k ( L(t 2 )—L(ti) ) ^> 1. The upper limit of the integral over k in Eq. [I] can then be cut off at 
2a7r/L(t 2 ) where a is a constant of 0{1) |TEf. For L(t 2 ) ^> L(ti), only the small k behavior 
of S(k,ti) contributes to the integral. Assume that linu^o S(k, t±) ~ k@ (/3 > 0). For 
quenches to zero temperature, S{k,t 2 ) will have the scaling form S(k,t 2 ) = L(t 2 ) d f(kL(t 2 )). 
Substituting into Eq. (P (with the appropriate limits of integration), gives 

i-2an/L(t 2 ) 



lim C(t u t 2 ) ~ L(t 2 y x < L(t 2 f 2 / dk k d ~ x k^ 2 f(kL(t 2 )), 

2>*1 JO 

^L{t 2 )-^' 2 . 



This immediately gives a lower bound on A, 

~ 2 
3 



The argument just presented, holds for conserved and nonconserved, scalar and vector order 
parameters. 

We now consider specific dynamical scenarios. Let Tj and Tp be the temperatures of 
the initial and final states respectively. We first focus on quenches from the high temper- 
ature phase (Tj = oo) to zero temperature (T F = 0). Since the initial state is disordered, 
linifc^o S(k, 0) ~ k°. In the absence of a conservation law, limfc_>o S(k, t\) ~ k° for both 
t\ = and ti in the scaling regime. Therefore (3 = and 

A > d/2. (2) 

This inequality was also obtained by Fisher and Huse using general scaling arguments 
and is consistent with all results to date P,p|,[7|-|lC|| . For ti = 0, conservation of the order 
parameter does not affect this inequality since f3 = for t± = 0. However, if t± is in the 
scaling regime, then lim^o S(k, t) ~ k 4 for d > 2 ]rj|] and (3 = 4. For d = 1, the dynamics 



is dominated by noise and Majumdar et al. find that linu^o S(k, t) ~ k 2 |TTJ, so that (3 = 2 
for d = 1. Therefore, for t\ in the scaling regime, 



| + 2 if d > 2, 
A> (3) 

| if d = l. 

These bounds suggest that the asymptotic exponent may depend on whether ti is, or is not 
in the scaling regime but do not rule out that the exponent is independent of t\. However 
for d = 1, Majumdar et al. find analytically and numerically that A = 1 for ti = 0, while we 
find that A > 3/2 for t\ in the scaling regime ]15| . 



For vector fields (with m, the number of components of the order parameter, > 2), an 
argument analogous to Ref. [pl|| , gives the same lim fc ^ S(k, t) ~ k A . This is supported by 



an extensive numerical integration of the Cahn-Hilliard equation | I6| . Therefore the lower 
bounds on A derived above are valid even for vector order parameters with m > 2. 

Quenches from the critical point (Tj = T c , Tp = 0) lead to long-range correlations of the 
initial configurations. In this case, A > d/2 no longer holds. More generally if S(k, 0) ~ k~ a 
we obtain A > (d — a)/2. (For critical dynamics o = 2 — rj where rj is the static critical 



exponent). This is consistent with the result of Bray et al. who found that, for nonconserved 
order parameter, A = (d — er)/2 for a greater than a critical value a c |J. 

Analysis of the bounds on the autocorrelation exponent for quenches to the critical point 
(T/ = oo, Tp = T c ), has to start afresh from Eq. (|I|). Since is in the critical point scaling 
regime, the correlation function has the following scaling form, S^/c, t 2 ) ~ k~ 2+r) f c (kL(t2)) . 
Substituting this form into Eq. (P gives A > (2d — 2 + rj + (3)/ 2. Therefore when t\ = we 
get A > (2d — 2 + T])/2. When t\ is also in the scaling regime, the bound on A depends on 
the behaviour of the scaling function f c (kL(t\)) as kL(t\) — > 0. For nonconserved systems 
lim x ^ fc(x) — > const. , or /3 = — 2 + 77 leading to A > d — 2 + r/. 

These lower bounds on A of course do not fix the value of the exponent. As previously 
mentioned, exact analytical and numerical computations on the 1-dimensional scalar model 
have been carried out for the case when t\ = 0. In higher dimensions, however, the empirical 



results are not very conclusive |L7 |. We therefore compute the asymptotic value of A by 



numerically integrating the Cahn-Hilliard equation in two- dimensions, 

^M = vVM), (4) 

where /i = —<fi + <fi 3 — V 2 0. We used an Euler discretization with St = 0.1 and 5x = 1.09 
and periodic boundary conditions. We discretize the Laplacian as 



V 2 , 



1 vg 
iij 6x 2 l + V2 



n.n.n. n.n. 



This choice decreases lattice anisotropy effects and allows a larger 5t before the onset of the 
checkerboard instability This dynamical equation is solved subject to random initial 



conditions which are uncorrelated and uniformly distributed between —0.05 and 0.05 (the 
initial state is disordered). Decreasing 5t has no effect on the numerical results. Increasing 
5x to 1.32 results in pinning effects which lead to a slower decay of the autocorrelation 
function at late times (even though the effect on the single-time behavior is less apparent). 

We have used the interfacial area density as a measure of the characteristic lengthscale 
L(t). Operationally, this is defined as (2 5xn x n y )/n opp , where n x n y is the total number of 



lattice sites and n opp is the number of sites with a nearest neighbor with of opposite sign. 
We recover the standard result that L(t) grows as t 1 / 3 for all t > 400. Other measures of 
the characteristic lengthscale, such as the first zero of the real space correlation function, 
also behave in the same manner (for t > 400). 

Fig. [I] shows C(t\,t) vs. L(t) for ti = and t between 100 and 12800 for three lattice 
sizes n = n x = n y = 64 (3084 initial conditions), n = 256 (1120 initial conditions) n = 
1024 (42 initial conditions). Concentrating on the largest lattice size, C(0,t) ~ L~ 3,7 for 
approximately two decades of t in its decay. There is a crossover to a slower decay at late 
times with C(0, t) ~ L~ 3 ' . At extremely late times there is an indication of an even slower 
decay, which we attribute to finite size effects. 

To emphasize the asymptotic trend, Fig. || shows L 3 C(0,t) vs. L(t) for the same data. 
Here it is clearer that the slower late time decay occurs at earlier times for smaller lattices, 
indicating finite size effects. The importance of the finite size effects was initially surprising 
since, for single time quantities, finite size effects only become important when L(t) is of 
order of the lattice dimension, L Q . Thus the usual length scales extracted from single-time 
quantities were identical for n = 256 and n = 1024. However, since C(0,t) decays rapidly 
with t, any small systematic effect becomes increasingly relevant as t increases. Clearly 
finite size effects on C(0,t) can be important (though not necessarily so) when the spread 
in C(0, t) is of the same order as C(0, t). The spread in C(0, t) decreases as L d ^ 2 and, 
based on our simulations, depends only weakly on L(t). Hence finite size effects can become 
important when C(0, t) ~ L(t)~ x ~ -^o^ 2 , i-e., much earlier than for single-time quantities. 
^From Fig. 0, C(0, t) for n = 64 first shows significant differences from the n = 256 result at 
C ~ 2.5 x 10~ 4 or L m 18. Thus we expect the n = 1024 data to be free of finite size effects 
down to C ~ 2.5 x 10 _4 /16 ~ 1.5 x 10~ 5 or L w 50. Therefore we make the preliminary 
conclusion that the true asymptotic value of A is approximately 3. Finally note that the 
late time result for n = 256 is consistent with X — d — 2. However, comparing this with the 
result for n = 1024 indicates that this regime is due to the finite size of the lattice. 

Fig. | shows C(t 1 ,t 2 ) vs. L(t 2 )/L{ti) for t\ = 100, 200 and 400. Although we cannot rule 
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out a further slower decay, our result is consistent with A ~ 4 for t\ in the scaling regime. 
Hence, we find that, for the d = 2, conserved, scalar model, A ~ 3 for t\ — and A ~ 4, for 
t\ in the scaling regime. To pin down these values more precisely would require simulations 
on larger lattices and for larger values of t\ and £2- 

Our numerical results for d = 2 and the lower bound on A for t\ in the scaling regime 
(Eq. (fl)) are inconsistent with the recent conjecture (Refs. |Tl|Jl2l ) that conservation of order 
parameter requires that A = d in all cases. We believe the apparent inconsistency is because 
the argument leading to A = d, incorrectly applies a scaling analysis to the non-scaling, 
k = mode. To be explicit, we briefly review the argument of Ref. [|TTJ (the argument 
presented in Ref. |L2| is similar). The Cahn-Hilliard equation (Eq. (^)) in k space is 

where D(k,t) = k 2 ^ k (t). Define S(k,0,t) = ( k (O)0_ k (t) ) = S(k,0,t) + V5 kfi (m 2 ). The 
formal solution to S(k, 0, t) is 

. >y(LV),kL(lf))\ 



S(k,0,t) = S(k, 0,0) exp (f dti 



h t> J 7 

where 7(L(t), kL{t)) = t ( D(k, t)0_k(t) ) /S(k, 0, t). In the scaling regime, 7(L, kL) = ^(kL) 
and we obtain (in the limit t ^> 0) , 



/' 

Jo 



l(kL(t')) f kL ^, 
Cll — 2/ I ctoc 

V JfcL(0) X 



zi(0) log (^y)+*7i (*£(*)) 



where L(t) ~ t l l z and 70 = lim 2 .^o+ 7( a; )- The result is 

S(k,0,t) = S(k,0,0)(^J F^kLit)). 

Since /dk S(k,0,t) = (0(r,O)0(r,t) ) ~ £(t)~ A = L(i)~ d+ ~< oz or A = d- 7o 2. However, since 
5*(/c = 0,0, t) is constant in time (conservation law !), it is argued that 70 must vanish and 
hence A = d. 

Our contention is that, even though the relation A = d — 70-2 with 70 = \\m x ^ + j(x) 
holds, the conclusion that, due to conservation of order parameter, 70 necessarily vanishes 



does not. This is because the relation X = d — ^qz is based on scaling, which only holds 
for k > 0, while the vanishing of 70 is based on the conservation of order parameter which 
only holds for k = 0. Therefore the conclusion 70 = is invalid since it applies a scaling 
argument to the nonscaling k = behavior. 

To see this more clearly, consider the quasi-static scattering intensity S(k,t) = S(k,t,t). 
For large t, S is given by S(k, t) = L(t) d f(kL(t)) + V5k,o{ m\ ). Global conservation of order 
parameter requires that S(k = 0,t) = constant, while the locally conservative dynamics 
requires that f(0) = 0. Hence global conservation leads to a discontinuity at k = for 
both S(k,ti,t 2 ) and j(L, kL) |fL9| . The singularities in 7(2) and S(k,ti,t 2 ) at k = can 
be removed by choosing the initial distribution so that ( m\ ) = 0. However, in this case, 
S(0, ti, t 2 ) vanishes independent of the value 0/70, so that the application of the conservation 
law does not fix the value of 70. (A. J. Bray has made a similar argument proving A = d 
does not necessarily hold p0|). 

Having provided useful lower bounds on A, we now ask whether it is possible to bound A 
from above? Unfortunately, we have not been able to provide useful upper bounds. However, 
we note that the bound (Eq. ([3])) as well as our numerical results violate the upper bound 
conjectured by Fisher and Huse, A < d 0. As they originally noted, this conjecture contains 
many assumptions. Moreover, in as much as their argument is aimed at the decay of the 
magnetization, their conjecture has validity only when the order parameter is not conserved 
and when ti = (so that S(k,ti) ~ A; ). 

We thank David Jasnow for the use of his workstations. We are grateful to John Ross 
and the University of Toronto Instructional and Research Computing Centre (UTIRC). This 
work was partially supported by the NSERC Canada. 
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FIGURES 

FIG. 1. The autocorrelation function C(0,t) for n = 64 (3084 initial conditions), 256 (1120 
initial conditions) and 1024 (42 initial configurations). Lines corresponding to (7(0, i) ~ L -3 and 
C(t) ~ L~ 2 are shown for comparison. 

FIG. 2. L 3 C(0,i) vs. L(t) for same data as Fig. 1. Lines corresponding to C(0, t) ~ L -3 and 
C(t) ~ L~ 2 are shown for comparison. These results emphasize that the slower decay at late times 
is due to finite size effects. 

FIG. 3. C(h, t 2 ) vs. L{t 2 )/L{ti) for n = 1024 and ti = 100, 200 and 400. Lines corresponding 
to C(ti,t2) ~ L~ 4 and C(ti,t2) ~ L~ 2 are shown for comparison. 
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We obtain bounds on the decay exponent A of the autocorrelation function in phase ordering dy- 
namics (defined by ]imt 2 >ti ( <f>(v, ti)<f>(r, t2)) ~ I(t2)~ X ). For non-conserved order parameter, we 
recover the Fisher and Huse inequality, A > d/2. If the order parameter is conserved we also find 
A > d/2 if t\ = 0. However, for t\ in the scaling regime, we obtain A > d/2 + 2 for d > 2 and A > 3/2 
for d = 1. For the one-dimensional scalar case, this, in conjunction with previous results, implies 
that A is different for t\ = and t\ ^ 1. In 2-dimensions, our extensive numerical simulations for 
a conserved scalar order parameter show that A fa 3 for t\ = and A fa 4 for t\ ^ 1. These results 
contradict a recent conjecture that conservation of order parameter requires A = d. Quenches to 
and from the critical point are also discussed. 



PACS: 64. 60. My, 68. 35. Fx 

Phase separation dynamics proceeds when a system is 
quenched from its high temperature, homogeneous phase 
to a low temperature, inhomogeneous phase (where sev- 
eral phases coexist in equilibrium). Due to its simple de- 
scription yet rich behavior, phase ordering dynamics has 
greatly enhanced our understanding of non-equilibrium 
processes [1]. At late times, the spatial distribution of 
domains can be described by a single time-dependent 
length, L(i) which typically grows algebraically in time, 
L(t) ~ t x l z . This results in a scale invariant equal- 
time correlation function C(r,t). More recently it has 
been realized that the unequal-time correlation function 
is scale covariant. In particular, the asymptotic decay 
of the two-time autocorrelation function, C(r, ti,t 2 ) = 
(<f)(r,ti)<f>(0,t2)} defines an independent exponent A, via 
lim tl<<2 C(0,ti,t 2 ) ~ (L(ti)/L(t 2 )) x ■ This exponent 
bears no apparent relation to the growth exponent z and 
so its value provides a sensitive test for approximate the- 
ories of phase ordering kinetics [2-8]. Although the auto- 
correlation function has been studied extensively for non- 
conserved order parameter dynamics [2,4-10], there has 
been hardly any work on conserved dynamics. However 
in a recent Letter, Majumdar et al. have shown numeri- 
cally and analytically that that A = 1 for m = 1, d = 1 
and t\ = where m is the number of components in the 
order parameter [11]. It has been further argued [11,12] 
that the conservation of order parameter demands that 
A = d for all m. 

In this Letter, we obtain lower bounds on the decay 
exponent A. For non-conserved order parameters, A > 
d/2 independent of t\, consistent with a general argument 
of Fisher and Huse [2]. For conserved order parameters, 
we also obtain A > d/2 for t\ = (assuming the quench 
is from a high temperature phase). However, for t\ in 
the scaling regime, we find that A > d/2 + 2 for d > 
2 and A > 3/2 for d = 1. This difference arises from 



the small k behavior of the scattering intensity S(k,ti). 
In conjunction with the exact result for A, for the 1- 
dimensional scalar model, with t\ =0 [11], we conclude 
that for d = 1, A depends on whether t\ = or t\ 1. 
To carry out the investigation in higher dimensions, we 
perform an extensive numerical integration of the Cahn- 
Hilliard equation (see Eq. (4) below) in d = 2. We find 
that A « 3 for t\ = and A « 4 for t\ in the scaling 
regime. This is inconsistent with the recent conjecture 
that A = d [11,12]. We discuss why this conjecture fails. 
We also derive bounds on A for quenches to and from the 
critical point. Our results easily extend to vector order 
parameters. 

We begin by obtaining the lower bounds on A. The 
equal point auto-correlation C(ti,t 2 ) = C(0,ti,t 2 ) is re- 
lated to the k space auto-correlation S^fc, ti , t 2 ) by 

C(*i,* 2 ) = y"dk(^k(*i)^-k(*2)> = f dk S(k,t 1 ,t 2 ). 

Here 4>(r,t) is the order parameter at point r and 
time t and 8(f)(r,t) = 4>(r,t) — mo with mo = 
V~ x ( dr <f)(r,t) and the fourier transform 6<f)k(t) = 
V -1 ' 2 J dr e~ lk r 6<f)(r,t). The angular brackets indicate 
an average over initial conditions. For a critical quench 
( mo) = and (m 2 ) is C(y _1 ), whereas for an off-critical 
quench ( mo ) is 0(1). 

Using the Cauchy-Schwartz inequality, we find 

c(t 1 ,t 2 )< j dk (<s<M*i)^-k(*i)) 1/2 <^2)^-k(* 2 ) 
~ J dk s , (Mi) 1/2 S'(M2) 1/2 . 

where S(k,t) = S(k,t,t). 

Now assume t 2 to be in the scaling regime with t 2 > ti ■ 
At late times, the scattering is due to the sharp inter- 
faces or defects. The k modes, 8(f)^ at times t\ and t 2 
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FIG. 1. The autocorrelation function C(0,t) for ra = 64 
(3084 initial conditions) 256 (1120 initial conditions) and 
1024 (42 initial configurations). Lines corresponding to 
C(0,t) ~ L~ 3 and C(t) ~ L~ 2 are shown for comparison. 



FIG. 2. £ 3 C(0,i) vs. L(t) for same data as Fig. 1. Lines 
corresponding to C(0, t) ~ i -3 and C(t) ~ i -2 are shown for 
comparison. These results emphasize that the slower decay 
at late times is due to finite size effects. 



will be uncorrelated when the interfaces move a distance 
greater than 2ir/k so that S{k,t\,t2) decreases rapidly 
for k ( L(t2) — L(t\) ) 1. The upper limit of the in- 
tegral over k in Eq. 1 can then be cut off at 2<nr / L(t2) 
where a is a constant of 0(1) [13]. For £(^2) L(t\), 
only the small k behavior of S{k,t\) contributes to the 
integral. Assume that lim^o S(k, t\) ~ k 13 (j3 > 0). For 
quenches to zero temperature, S{k,t2) will have the scal- 
ing form S(k,t2) = L(t2) d f(kL(t2))- Substituting into 
Eq. (1) (with the appropriate limits of integration), gives 

lim C(t 1 ,t 2 ) ~L{t 2 y X 

[■2air I L(t 2 ) 

< L(t 2 ) d/2 / dk k*- 1 k p l 2 f[kHt 2 )).. 

Jo 

~ L(t 2 )- (d+ P )/2 . 
This immediately gives a lower bound on A, 



A > 



The argument just presented, holds for conserved and 
nonconserved, scalar and vector order parameters. 

We now consider specific dynamical scenarios. Let 
Tj and Tp be the temperatures of the initial and fi- 
nal states respectively. We first focus on quenches from 
the high temperature phase (Tj = 00) to zero temper- 
ature (Tp = 0). Since the initial state is disordered, 
limfc^o S(k, 0) ~ k° . In the absence of a conservation 
law, limfc^o S(k, ti) ~ k° for both t\ = and t\ in the 
scaling regime. Therefore [3 = and 



A > d/2. 



(2) 



This inequality was also obtained by Fisher and Huse us- 
ing general scaling arguments [2] and is consistent with 
all results to date [4,5,7-10]. For t\ = 0, conservation of 
the order parameter does not affect this inequality since 



f3 = for t\ = 0. However, if t\ is in the scaling regime, 
then lim^o S^M) ~ k 4 for d > 2 [14] and j3 = 4. For 
d = 1, the dynamics is dominated by noise and Majum- 
dar et al. find that lim^o S(k, t) ~ k 2 [11], so that (3 = 2 
for d = 1. Therefore, for t\ in the scaling regime, 



A > 



if, 
if, 



>2, 
= 1. 



(3) 



These bounds suggest that the asymptotic exponent may 
depend on whether t\ is, or is not in the scaling regime 
but do not rule out that the exponent is independent of 
t\. However for d = 1, Majumdar et al. find analytically 
and numerically that A = 1 for t\ = 0, while we find that 
A > 3/2 for t\ in the scaling regime [15]. 

For vector fields (with m, the number of components 
of the order parameter, > 2), an argument analogous to 
Ref. [14], gives the same lim^o S(k, t) ~ k 4 . This is 
supported by an extensive numerical integration of the 
Cahn-Hilliard equation [16]. Therefore the lower bounds 
on A derived above are valid even for vector order param- 
eters with m > 2. 

Quenches from the critical point (T/ = T c ,Tp = 0) 
lead to long-range correlations of the initial configura- 
tions. In this case, A > d/2 no longer holds. More gen- 
erally if S(k,0) ~ k~ a we obtain A > (d - a)/2. (For 
critical dynamics a = 2 — rj where rj is the static critical 
exponent). This is consistent with the result of Bray et 
al. who found that, for nonconserved order parameter, 
A = (d — <j)/2 for it greater than a critical value <r c [6]. 

Analysis of the bounds on the autocorrelation expo- 
nent for quenches to the critical point (T/ = 00, Tp = T c ), 
has to start afresh from Eq. (1). Since ^2 is in the criti- 
cal point scaling regime, the correlation function has the 
following scaling form, S(k,t2) ~ k~ 2+v f c (kL(t2)). Sub- 
stituting this form into Eq. (1) gives A > (2d— 2+r)+/3)/2. 
Therefore when t\ = we get A > (2c? — 2 + rj)/2. When 
t\ is also in the scaling regime, the bound on A depends 



2 



on the behaviour of the scaling function f c (kL(ti)) as 
kL(t\) — ► 0. For nonconserved systems lirm^o fc(%) —> 
const. , or [3 = —2 + r\ leading to A > d — 2 + rj. 

These lower bounds on A of course do not fix the value 
of the exponent. As previously mentioned, exact ana- 
lytical and numerical computations on the 1-dimensional 
scalar model have been carried out for the case when 
ti = 0. In higher dimensions, however, the empirical re- 
sults are not very conclusive [17]. We therefore compute 
the asymptotic value of A by numerically integrating the 
Cahn-Hilliard equation in two-dimensions, 



d<j>(r,t) 
dt 



V 2 /i(r,i), 



(4) 



■f <f> 3 — V 2 </>. We used an Euler discretiza- 
0.1 and Sx = 1.09 and periodic boundary 



where \i = - 
tion with St 
conditions. We discretize the Laplacian as 



1 



V2 



1 + V2 



This choice decreases lattice anisotropy effects and al- 
lows a larger St before the onset of the checkerboard in- 
stability [18]. This dynamical equation is solved subject 
to random initial conditions which are uncorrelated and 
uniformly distributed between —0.05 and 0.05 (the ini- 
tial state is disordered). Decreasing St has no effect on 
the numerical results. Increasing Sx to 1.32 results in 
pinning effects which lead to a slower decay of the auto- 
correlation function at late times (even though the effect 
on the single-time behavior is less apparent). 

We have used the interfacial area density as a measure 
of the characteristic lengthscale L(t). Operationally, this 
is defined as (2 Sx n x n y )/n opp , where n x n y is the total 
number of lattice sites and n opp is the number of sites 
with a nearest neighbor with <f> of opposite sign. We re- 
cover the standard result that Lit) grows as t 1 ! 3 for all 
t > 400. Other measures of the characteristic length- 
scale, such as the first zero of the real space correlation 
function, also behave in the same manner (for t > 400). 

Fig. 1 shows C(ti,t) vs. L(t) for t\ = and t between 
100 and 12800 for three lattice sizes n = n x = n y = 64 
(3084 initial conditions), n = 256 (1120 initial condi- 
tions) n = 1024 (42 initial conditions). Concentrating on 
the largest lattice size, C(0,t) ~ L~ 37 for approximately 
two decades of t in its decay. There is a crossover to a 
slower decay at late times with C(0,t) ~ L~ 30 . At ex- 
tremely late times there is an indication of an even slower 
decay, which we attribute to finite size effects. 

To emphasize the asymptotic trend, Fig. 2 shows 
L 3 C(0,t) vs. L{t) for the same data. Here it is clearer 
that the slower late time decay occurs at earlier times 
for smaller lattices, indicating finite size effects. The im- 
portance of the finite size effects was initially surprising 
since, for single time quantities, finite size effects only 
become important when L(t) is of order of the lattice 
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log 10 (L(t 2 )/L( tl )) 

FIG. 3. C(t 1 ,t 2 ) vs. £(i 2 )/£(ii)for n = 1024 and ii = 100, 
200 and 400. Lines corresponding to C(t\,t2) ~ £ -4 and 
C(t\,t2) ~ L~ 2 are shown for comparison. 



dimension, Lq. Thus the usual length scales extracted 
from single-time quantities were identical for n = 256 and 
n = 1024. However, since C(0,t) decays rapidly with t, 
any small systematic effect becomes increasingly relevant 
as t increases. Clearly finite size effects on C(0,t) can be 
important (though not necessarily so) when the spread 
in C(0,t) is of the same order as C(0,t). The spread in 
C(0,t) decreases as L d ^ 2 and, based on our simulations, 
depends only weakly on L(t). Hence finite size effects can 
become important when C(0,t) ~ L(t)~ x ~ L Q d ^ 2 , i.e., 
much earlier than for single-time quantities. ^From Fig. 
2, C(0,t) for n = 64 first shows significant differences 
from the n = 256 result at C ~ 2.5 x 10~ 4 or L k, 18. 
Thus we expect the n = 1024 data to be free of finite 
size effects down to C ~ 2.5 x 10 _4 /16 « 1.5 x 10~ 5 or 
L « 50. Therefore we make the preliminary conclusion 
that the true asymptotic value of A is approximately 3. 
Finally note that the late time result for n = 256 is con- 
sistent with A = d = 2. However, comparing this with 
the result for n = 1024 indicates that this regime is due 
to the finite size of the lattice. 

Fig. 3 shows C(t 1; t 2 ) vs. L(t 2 )/L(t 1 ) for t 1 = 100, 
200 and 400. Although we cannot rule out a further 
slower decay, our result is consistent with A « 4 for t\ in 
the scaling regime. Hence, we find that, for the d = 2, 
conserved, scalar model, A « 3 for t\ = and A « 4, for 
t\ in the scaling regime. To pin down these values more 
precisely would require simulations on larger lattices and 
for larger values of t\ and t 2 - 

Our numerical results for d = 2 and the lower bound 
on A for t\ in the scaling regime (Eq. (3)) are inconsistent 
with the recent conjecture (Refs. [11,12]) that conserva- 
tion of order parameter requires that A = d in all cases. 
We believe the apparent inconsistency is because the ar- 
gument leading to A = d, incorrectly applies a scaling 
analysis to the non-scaling, k = mode. To be explicit, 
we briefly review the argument of Ref. [11] (the argu- 
ment presented in Ref. [12] is similar). The Cahn-Hilliard 



3 



equation (Eq. (4)) in k space is 



dt 



D(k,t), 



where D(k,t) = k 2 ^ k {t). Define S{k,Q,t) = 
(^k(O)^-k(O) = S(k,0,t) + V8 kt o(m 2 ). The formal so- 
lution to S(k, 0, t) is 



S(Jfe,(M) = S(M,0)exp( /* dt' iWDlM^l 



where j(L(t), kL(t)) = t ( D(k, t)<f>_ k (t) )/S(k, 0, t). In 
the scaling regime, j(L,kL) = j(kL) and we obtain (in 
the limit t > 0) , 



L(0) 



27(0) log 



z 7 i (*£(*)) : 



where ~ t 1 ^ and 70 = lim a ._ i . + tM- The result is 
S(k,0,t) = S(k,0,0) (|||) 70 F.ikLit)). 

Since Jrfk 5(A,0,<) = ( <f>(r, 0)<f>(r, t) } ~ £(*)~ A = 
or A = d - j z. However, since S(k = 0, 0, t) 
is constant in time (conservation law !), it is argued that 
70 must vanish and hence A = d. 

Our contention is that, even though the relation A = 
d — Joz with 70 = linx,,^o+ l{ x ) holds, the conclusion 
that, due to conservation of order parameter, 70 neces- 
sarily vanishes does not. This is because the relation 
A = d — 70Z is based on scaling, which only holds for 
k > 0, while the vanishing of 70 is based on the conser- 
vation of order parameter which only holds for k = 0. 
Therefore the conclusion 70 = is invalid since it applies 
a scaling argument to the nonscaling k = behavior. 

To see this more clearly, consider the quasi-static scat- 
tering intensity S{k,i) = S{k,t,i). For large t, S is given 
by S(k,t) = L(t) d f(kL(t)) + Vhfi{ml). Global conser- 
vation of order parameter requires that S(k = 0, t) = con- 
stant, while the locally conservative dynamics requires 
that /(0) = 0. Hence global conservation leads to a dis- 
continuity at k = for both S(k, t\, t'2) and j(L, kL) [19]. 
The singularities in j(x) and S(k,ti,t2) at k = can 
be removed by choosing the initial distribution so that 
(mg) = 0. However, in this case, S(0,ti,t2) vanishes 
independent of the value 0/70, so that the application of 
the conservation law does not fix the value of 70. (A. J. 
Bray has made a similar argument proving A = d does 
not necessarily hold [20]). 

Having provided useful lower bounds on A, we now ask 
whether it is possible to bound A from above? Unfor- 
tunately, we have not been able to provide useful upper 



bounds. However, we note that the bound (Eq. (3)) as 
well as our numerical results violate the upper bound 
conjectured by Fisher and Huse, A < d [2]. As they orig- 
inally noted, this conjecture contains many assumptions. 
Moreover, in as much as their argument is aimed at the 
decay of the magnetization, their conjecture has valid- 
ity only when the order parameter is not conserved and 
when ti = (so that S(Mi) ~ 
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